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ABSTRACT 


The  mathematical  model  presently  used  by  the  Space  Detection  and 
Tracking  System  (SPADATS) ,  operated  by  the  U.S.  Air  Force,  Air  Defense 
Command,  to  describe  the  orbital  motion  of  Earth  satellites  is  presented. 

The  formulation  is  given  and  the  non-Keplerian  terms  are  justified. 

Various  output  formats  of  the  system  are  described  and  defined.  Limitations 
of  the  model  are  discussed  in  a  qualitative  fashion. 
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THE  SPADATS  MATHEMATICAL  MODEL 


"Glendower:  I  can  call  spirits  from  the  vasty  deep. 

Hotspur:  Why,  so  can  I  and  so  can  any  nan; 

But  will  they  come  when  you  do  call  for 
then? " 

Henry  IV.  Part  1,  Act  3,  Scene  1 


Anybody  can  predict  the  future  positions  of  Earth  satellites. 
Whether  they  will  be  at  the  predicted  place  is  another  natter.  In  nore 
scientific  language:  the  problem  is  one  of  accuracy  in  prediction.  This 
accuracy  or,  more  realistically,  the  inaccuracy  is  due  to  many  causes. 

The  purpose  of  this  paper  is  to  discuss  various  mathematical  models  of 
satellite  orbits,  the  underlying  assumption,  the  errors  and  limitations 
of  the  models.  In  particular,  the  SPADATS  model  will  be  defined  and 
delimited. 


THE  NATURE  OF  MATHEMATICAL  MODELS 

A  mathematical  model  is  an  entity  in  itself.  It  does  not 
need  any  physical  theory  to  sustain  it.  The  test  of  the  model  is 
whether  it  generates  predictions  of  the  required  accuracy.  For  instance, 
by  a  few  days  of  observation  of  the  rising  Sun,  it  is  possible  to  predict 
where  on  the  horizon  it  will  rise  the  next  day.  After  records  have  been 
kept  for  about  a  y<ur,  it  should  be  possible  to  design  a  mathematical 
model  for  the  azimuth  of  the  rising  Sun  (at  a  particular  site)  which 
Includes  an  annual  term.  This  model  would  be  quite  accurate  for  many 
years ,  if  we  know  how  long  a  year  is . 

This  is  an  important  point.  The  accuracy  of  the  mathematical 
model  depends  on  the  accuracy  of  its  constants.  It  may  take  many  years  of 
counting  the  days  between  the  risings  of  the  Sun  at  a  given  place  on  the 
horizon  to  refine  the  value  to  365  days  and  the  best  part  of  a  century  to 
refine  it  to  365.24.  Why?  Because  the  measurement  of  azimuth  will  be 
inaccurate.  It  is  necessary  to  define  the  given  place  on  the  horizon 
accurately  and  the  place  to  stand.  Thus,  the  accuracy  of  the  constants 


depends  on  Che  accuracy  of  the  measurements  used  Co  obtain  the  constants. 

Such  a  mathematical  model  was  arrived  at,  many  millenia  ago. 
without  the  knowledge  of  any  underlying  causes.  It  is  possible  to  build 
quite  an  elaborate  structure  this  way. 

The  Ptolemaic  System  of  Planetary  Motions  was  based  on  rather 
unsound  principles.  Yet,  the  constants  wore  %*?11  determined  and  the 
system  represented  the  phenomena  quite  accurately.  The  more  realistic 
universe  of  Aristarchos  (fl.  280  B  C.)  with  the  Sun  at  its  center  was 
rejected  partly  because  Aristarchos  had  to  place  the  stars  almost  infinitely 
far  away,  which  seemed  ridiculous,  and  partly  because  hardly  any  simplifi¬ 
cation  could  be  achieved  until  the  Pythagorean  assumption  of  circular 
motions  was  overthrown  by  Kepler. 


THE  EARTH -SATELLITE  PROBLEM 

The  work  of  Newton  makes  it  possible  to  express  the  problem  of 
the  prediction  of  the  future  positions  of  a  body  in  terms  of  the  forces 
acting  on  it.  One  may,  therefore,  integrate  the  resultant  accelerations 
numerically  or  analytically  to  obtain  positions  and  velocities. 

The  dominant  force  on  Earth  satellites  (except  during  launch  or 
re-entry)  is  the  central  attraction  of  the  Earth.  If  this  were  the  only 
force,  the  orbit  of  the  satellite  would  be  a  Keplerian  ellipse.  This  type 
of  motion  is  completely  defined  by  the  six  elements  which  are  the  constants 
of  an  analytical  integration. 

The  elements  are  mathematically  equivalent  to  the  position  and 
velocity  of  the  vehicle  at  any  given  time.  Three  elements  describe  the 
shape  and  size  of  the  orbit  and  the  location  of  the  vehicle  in  the  orbit 
at  a  specified  time  The  classical  dimensional  elements  are  as  follows: 


a  -  semi-major  axis 
e  -  eccentricity 

T  -  time  of  perigee  passage  (i.e.,  point  of  closest  approach 
to  the  Earth) 


The  remaining  three  elements  specify  the  orientation  of  the  conic 
in  space.  The  most  descriptive  elements  are  the  orientation  angles.  These 
angles  are  as  follows: 

1  -  inclination  of  the  orbit  with  respect  to  the  equa~ 

torial  plane 

Cl  -  right  ascension  of  the  ascending  node,  the  angle 

in  the  equatorial  plane  between  the  vernal  equinox 
and  the  line  to  the  ascending  node 

U)  -  argument  of  perigee,  the  angle  in  the  orbit  plane 

between  the  line  to  the  ascending  node  and  the  line 
through  perigee . 


The  orientation  angles  and  some  related  unit  vectors  are  shown  in  Figure  1. 
In  actual  computation,  it  is  preferable  to  use  the  components  of  those 
vectors  (direction  cosines)  directly  and  to  disregard  the  angles  as  much 
as  possible. 

There  are,  however,  additional  forces  acting  on  the  satellite. 
These  may  be  treated  as  perturbations  on  the  Kepler  solution.  In  this 
sense  it  may  be  said  that  the  orbit  of  an  Earth  satellite  is  continually 
changing  due  to  perturbative  effects.  Though  much  smaller  than  the  central 
attraction  of  the  Earth,  the  perturbative  forces  can  cause  deviations  from 
the  basic  Keplerian  ellipse  which  are  appreciable  from  the  standpoint  of 
satellite  tracking.  Some  of  the  causes  that  may  be  important  in  terms  of 
the  SPADATS  mission  are: 


(1)  Earth  bulge 

(2)  Atmospheric  drag 

(3)  Solar  attraction  and  radiation  pressure 
(A)  Lunar  gravity 


The  effects  of  the  perturbations  are  usually  grouped  into  two 
classes;  secular  (those  which  grow  with  time)  and  periodic  (those  which 
oscillate).  The  secular  class  causes  the  larger  errors,  because  these 
perturbations  build  up  with  time.  In  some  cases  and  applications,  the 
periodic  perturbations  of  an  orbit  cannot  be  ignored. 
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To  fulfill  Che  mission  of  satellite  tracking,  the  future  courses 
of  satellites  must  be  predicted  with  some  degree  of  precision.  The  pertur¬ 
bative  accelerations  may  be  either  integrated  numerically  or  by  the  term-by 
term  analytical  integration  of  a  series  expansion.  The  former  process  is 
called  Special  Perturbations  in  Celestial  Mechanics.  The  analytic  integra¬ 
tion  results  in  a  General  Perturbations  •theory*’  of  the  orbit.  One  advan¬ 
tage  of  such  a  "theory*1  involves  the  ability  to  predict  the  position  for  any 
time  without  having  to  compute  intermediate  positions  as  is  necessary  for 
numerical  integration.  Again,  the  gross  effects  of  the  perturbative  force 
become  more  readily  apparent  from  the  theory  than  from  inspection  of  the 
numbers  resulting  from  Special  Perturbations  calculations.  In  the  Special 
Perturbations  process  the  perturbative  expressions  are  generally  simpler  to 
express  and  formulate  with  comparable,  or  better,  accuracy. 

The  ultimate  in  simplicity  results  when  even  the  conic  solution 
is  ignored  and  the  total  acceleration  is  integrated  numerically.  This  pro¬ 
cess  is  called  "Cowell’s  Method"  to  distinguish  it  from  the  older  methods 
l.e.,  Variation  of  Elements  and  Encke's  Method.  Cowell's  method  became 
feasible  with  the  development  of  desk  calculators  and  is  especially  suit¬ 
able  for  large-scale  electronic  computers;  however,  it  is  relatively  slow 
for  satellite  orbits  unless  the  perturbative  acceleration  is  almost  equal 
to  the  gravitational  acceleration. 

General  Perturbations  theories  have  been  developed  to  express  the 
effect  of  the  aspherical  gravitational  field  of  the  Earth,  which  causes  the 
principal  perturbations  on  a  near-Earth  satellite  during  roost  of  its  orbital 
lifetime.  The  SPADATS  mathematical  model  is  based  mainly  on  the  work  of 
Kozai  and  Brouwer  (See  Appendix  C) . 

It  is  possible  to  make  an  increasingly  accurate  Integration  of  the 
equations  of  motion  in  specific  cases;  however,  it  is  usually  unwarranted  to 
complicate  the  solution  because  the  apparent  satellite  path  will  deviate 
from  the  true  path  locauae  of  inexact  observations,  inexact  theory  and 
constants  (model  errors),  and  inexact  computations  (e.g.,  round-off).  The 
majority  of  the  Earth  satellite  prediction  needs  can  be  met  practically  by 
a  relatively  simple  General  Perturbations  treatment  in  conjunction  with  a 
correction  procedure  which  employs  observational  data.  This  correction  is 
an  integral  part  of  any  computation  of  an  actual  orbit  because  it  establishes 
the  elements  of  the  orbit.  The  form  and  size  of  the  error  in  the  elements 
depends  not  only  on  the  observations  used  in  the  correction  but  somewhat  on 
the  correction  process  also. 
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THE  DIFFERENTIAL  CORRECTION 


In  order  to  correct  an  orbit,  it  is  necessary  to  '‘represent’-  the 
observations  (i.e.,  calculate  what  the  observations  would  have  been  had  the 
elements  and  the  mathematical  model  for  the  orbit  representation  been  cor¬ 
rect).  The  end  products  of  the  observation  processing  are  the  residuals 
which  are  the  differences:  the  observed  quantities  minus  the  'represented' 
quantities . 

The  calculation  of  the  represented  observation  may  include  all 
components  of  a  complete  'fir.''  (three  components  of  position  and  three 
components  of  velocity).  However,  residuals  can  be  obtained  for  only  those 
components  which  are  actually  observed.  Some  of  these  residuals  may  not 
be  usable  if  the  observing  instrument  performs  poorly  in  some  aspects  of 
measuring.  Nevertheless,  the  remaining  residuals  may  form  the  basis  for 
differentially  correcting  an  orbit. 

Uhen  there  is  redundancy  in  the  observations,  i.e.,  when  there  are 
more  residuals  than  there  are  parameters  to  be  corrected  (normally  seven), 
the  differential  correction  process  becomes  even  more  effective,  since  the 
process  makes  use  of  the  method  of  statistical  averaging  known  as  least 
squares.  Observations  which  are  separated  by  long  intervals  of  time  and 
which  are  obtained  from  different  points  on  the  Earth  can  be  included. 

The  observations  are  weighted,  in  the  present  system,  according 
to  certain  geometric  effects.  These  weighting  factors  are  given  in  Table  I. 
Except  for  the  range-rate  residual,  these  factors  express  the  physical 
displacement  (in  kilometers)  of  the  vehicle  from  its  assumed  place  as 
indicated  by  the  observed  quantity.  This  physical  displacement  component 
is  also  used  to  reject  those  equations  of  condition  which  indicate  a  dis¬ 
placement  greater  than  an  absolute  distance  criterion  or  greater  than 
1.5  times  the  root-means -square  (r.m.s.)  of  the  displacement  component.* 

The  equations  for  lange-rate  are  tested  against  a  separate  but  similar 
absolute  criterion  and  against  1.5  times  the  r.m.s.  of  the  remaining 
range-rate  residuals.  This  rejection  process  may  be  considered  a  0  or  1 
weighting  on  the  basis  of  the  internal  consistency  of  the  data.  The  factors 
cos  h  and  cos  6  account  for  the  degradation  of  azimuth  and  right  ascension 
measures  near  the  zenith  and  the  poles,  respectively.  Weighting  on  the 
basis  of  sensor  characteristics  is  not  performed  at  present. 


♦remaining  after  rejection  or  the  basis  of  the  absolute  criterion. 
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TABLE  I 


MULTIPLICATIVE  FACTORS  OF  EQUATIONS  OF  CONDITION 


Observed  Quantity 


Factors 


Residual  Form 


Range 
Range-Rate 
Azimuth 
Elevation 
Right  Ascension 


None  Ap 

Range  PAp 

Range,  Cosine  Elevation  Pc os  h  A  a 

Range  pA  h 

Range,  Cosine  Declination  Pcos  6Aa 

pA6 


Declination 


Range 


THE  SPADAT  SYSTEM  ORBIT  REPRESENTATION 


The  angles  which  are  used  in  the  classical  element  set  have  the 
disadvantages  of  singularities  and  computational  inef f iencies .  For  example, 
when  the  eccentricity  becomes  small,  the  argument  andQtii9e  of  perigee  become 
ill-defined.  When  the  inclination  is  small  (nearly  0  or  nearly  180°),  the 
node  is  ill-defined.  As  usual,  the  exact  singularities  have  orbit  represen¬ 
tations  which  are  often  simpler  than  those  for  the  general  conic.  However, 
the  near-singularities  may  cuase  difficulties  in  both  representation  and 
correction. 

Some  of  these  difficulties  can  be  avoided,  and  computational 
efficiency  can  be  gained,  by  use  of  an  element  set  based  on  vectors  (direc¬ 
tion  cosines)  Instead  of  angles.  Such  a  set  is  the  Element  Set  listed  below. 
The  relationships  between  the  orientation  vectors  and  the  orientation  angles 
are  illustrated  in  Figure  1. 


TABLE  II 

INTERNAL  ELEMENTS  OF  SPADAT  SYSTEM 


t  ■  epoch  time;  can  be  selected  at  an  arbitrary  position  in  the 

°  orbit 

axN,av n“  components  of  a  along  N  and  M,  respectively;  where  a  ■  eP, 
y  "  -  e  cos  U) ,  a  „  -  e  sin  U) 


h 

L 


xN 


yN 


P  is  a  unit  vector  directed  toward  perigee  from  the  center 
of  the  Earth 

^W,  angular  momentum  vector,  where  W  is  the  normal  to  the 
orbit  plane  (See  Figure  1)  and  p  is  the  semi-latus  rectum 

cos  i _ 


-  «  +  M  +  ,  . 

u  0  O  I  COS  1  I  o 

■  mean  longitude  at  epoch 

where  Mq, ,  fl  ,  are  the  mean  anomaly,  argument  of  perigee,  and  right  ascen¬ 
sion  of  the  node,  all  at  epoch 


c  ■  period  decay  rate 

d  -  period  decay  acceleration 


drag  parameters 
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The  formulation  presently  used  to  obtain  position  and  velocity 
from  these  elements  is  contained  in  Appendix  A.  The  model  includes  the 
first-order  secular  terms  in  L,to  and  ft,  due  to  the  second  zonal  harmonic 
of  the  geopotential,  as  well  as  the  large  long-period  variations  in  L  and 
a  due  to  the  third  zonal  harmonic.  The  derivation  of  the  formulas  in¬ 
volving  those  harmonics  is  given  in  Appendix  C.  The  definition  of  output 
elements  in  terms  of  those  given  in  Table  11  is  discussed  in  Appendix  B. 


LIMITATIONS 

The  principal  limitation  of  the  SPADATS  Model  is  in  accuracy.* 

To  preserve  computational  efficiency  only  a  limited  number  of  perturbative 
terms  are  included  in  the  representation.  No  short -period  effects  are  in¬ 
cluded.  The  elements  determined  by  the  differential  correction  are  mean 
elements  at  the  epoch,  t  ,  that  is,  they  are  intended  to  be  the  elements 
which  would  represent  the  orbit  if  all  the  periodic  terms  were  zero  at  tQ. 

The  elements  are  produced  by  the  SPADAT  center  in  various  forms. 
Several  of  these  forms  are  defined  in  Appendix  B.  One  form,  the  Space  Track 
Bulletin,  is  no  longer  widely  circulated.  It  consists  of  four  parts.  Part  I 
contains  some  of  the.  same  information  as  on  the  7-Element  Cards.  Part  II 
gives  the  times,  longitudes  and  heights  of  a  set  of  ascending  nodal  cross¬ 
ings.  Part  III  gives  the  time  and  longitude  corrections  to  convert  the 
information  from  Part  II  to  other  latitude  crossings.  Part  III  is  calculated 
for  some  median  revolution  of  Part  II.  Part  IV  is  the  SATOR  code,  giving  the 
elements  in  a  telegraphic  format.  Since  it  receives  wider  distribution,  the 
SATOR  code  is  also  defined  in  Appendix  B. 

It  must  be  borne  in  mind  that  the  elements  (and  the  bulletins)  are 
issued  for  purposes  of  prediction.  The  elements  are  determined  from  observa¬ 
tions  made  prior  to  -neir  epoch.  The  bulletin  elements,  for  example,  are 
best  for  that  preceding  period  and  not  for  the  revolutions  indicated  in  Partll. 
Furthermore,  they  are  the  mean  elements  from  a  truncated  model  and  cannot  be 
used  to  define  starting  conditions  for  a  numerical  integration  of  the  orbit. 


*Quantitative  information  on  the  accuracy  of  the  present  SPADAT  System  may 
be  found  in  Aeronutronic  Report  S-2071,  ’'Spadat  System  Observational  Data 
Studies,1'  Volume  I  of  Studies  of  Satellite  Prediction  Accuracy.  ESD-TDR-63-123 , 
by  G.  R.  Miczaika  and  D.  M.  Slager. 


-9 


It  may  be  thought  that  it  is  possible  to  obtain  osculating  elements 
by  adding  the  periodic  terms  to  the  mean  elements.  Such  osculating  elements 
may  also  be  inaccurate.  The  differential  correction  obtains  the  best  set  of 
mean  elements  from  the  observational  material  presented  to  it.  The  associated 
least-squares  process  is  designed  to  minimize  the  transfer  of  random  errors 
to  the  elements.  It  is  limited,  however,  by  the  absence  of  some  of  the  peri¬ 
odic  terms  in  the  representation.  Least-squares  would  be  able  to  smooth  out 
these  small  variations  if  the  observations  presented  it  with  a  sample  of 
these  effects  that  is  random  with  respect  to  phase.  The  sample  of  the 
short-periodic  effects  is  biased  by  the  preponderance  of  sensors  in  the  North¬ 
ern  Hemisphere  and  by  a  frequent  concentration  of  observations  at  the  northern 
apex  of  the  orbits.  The  sample  of  the  long-period  effects  is  biased  because 
of  the  limited  time-span  of  observations  used  in  a  typical  differential  correc¬ 
tion.  This  time-span,  dictated  usually  by  the  irregular  variations  of  atmos¬ 
pheric  drag,  may  not  allow  a  significant  arc  of  long-period  terms  to  be 
included. 


Assume  that  the  observations  are  so  concentrated,  for  one  of  the 
foregoing  reasons,  that  they  sample  some  neglected  periodic  term  in  a  small 
part  of  a  full  cycle  (See  Figure  2)  .  Then  the  SPADATS  "mean'*  element  or 
position  would  tend  to  include  the  mean  value  of  the  neglected  term  during 
the  observation  span.  The  figure  indicates  that  the  epoch  is  taken  at  the 
end  of  the  observation  span.  This  would  be  the  case  if  the  neglected  term 
had  a  long  period,  since  the  SPADAT  System  places  the  epoch  at  a  nodal  cross¬ 
ing  near  the  last  observation.  It  can  be  seen  that  adding  the  value  of  the 
neglected  term  to  the  SPADATS  mean  does  not  necessarily  give  the  osculating 
value  at  the  epoch  and  may,  Indeed,  give  a  more  inexact  value.  This  indirect 
effect  of  the  neglected  terms  is  somewhat  mitigated  by  including  in  the  dif¬ 
ferential  correction  different  types  of  observations  (e.g.,  range-rate  as 
well  as  range)  which  simultaneously  sample  the  same  periodic  effects  at  dif¬ 
ferent  phases . 

The  deviation  of  the  SPADATS  Model  from  the  true  orbit  is  kept  in 
bounds  by  testing  the  new  observations  (after  epoch) .  If  the  root -mean- 
square  of  the  displacement  of  the  observed  positions  from  the  positions  com¬ 
puted  from  the  model  exceeds  a  bound,  a  new  correction  is  called  for  by  the 
system. 


Another  type  of  limitation  concerns  the  applicability  of  the  semi¬ 
automatic  programming  system,  which  uses  the  appended  formulation.  This  sys¬ 
tem  is  only  designed  to  handle  near-Earth  satellites.  Separate  programs  must 
be  used  to  predict  the  position  of  lunar  probes,  space  probes,  and  equatorial 
satellites . 
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FIGURE  2.  EFFECT  OF  NEGLECTED  PERIODIC  TERM 


SAT0R  CODE 


The  SATOR  Code  is  produced  as  Part  IV  of  the  Bulletin  by  the  BLTNSGP 
program.  The  format  is  as  follows: 

SAT0R  aabbc  deeff  ggggZ  hhhhX  K0WES 

iiiii  jkkkk  ARPER  11111  mnnnX  PERI0D 

OOOOO  PPPPP  ECCEN  qqqqq  PERRA  rrrrr 

RAF  RE  sssss  (sssss  repeated  RTASC  ttttt 

as  necessary) 

Key: 

aa  ■  Last  two  digits  of  year  satellite  launched 
bb  ■  Greek  letter  designation,  01  -  Alpha, 

02  -  Beta,  etc.,  or  number  after  1963 
c  “  Component  A-l,  B«2,  etc. 
d  -  Reference  time  (time  of  perigee  passage 

closest  to  epoch) :  last  digit  of  numerical 
notation  for  month;  i.e.,  1* January  or 
November,  2«February  or  December,  3-March, 
etc. 

ee  -  Reference  time:  date 
ff  -  Reference  time:  hour 
gggg  -  Reference  time:  minutes  and  hundredths  of 
minutes 

Z  -  Universal  time,  Greenwich  Mean  Time 
hhhh  -  Inclination  in  degrees  and  hundredths  of 
degrees.  If  the  orbit  inclination  is 
negative  (satellite  fired  westward)  group 
is  preceded  by  NEGAT. 

X  ■  Always  an  X 

N0WES  -  Subindicator  for  geographical  longitude 
of  orthbound  node  west  of  Greenwich  at 
reference  time 

iiiii  ■  West  longitude  of  northbound  node,Xw>in 
degrees  and  hundredths  of  degrees, 
j  ■  1.  If  plus:  When  the  "prime  sweep  interval" 
is  one  day  plus  a  certain  number 
of  minutes. 

2.  If  minus: When  the  "prime  sweep  .interval" 

is  one  day  minus  a  certain  number 
of  minutes.  This  is  equivalent 
to  saying  that  the  same  portion 
of  the  orbit  plane  will  reappear 
at  the  same  location  a  certain 
number  of  minutes  earlier  each 
day . 
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kkkk  *  Number  of  minutes  by  which  "prime  sweep  interval" 
differs  from  one  day  or  1440  minutes.  This  is 
another  way  of  expressing  the  relative  motion  of 
the  orbit  plane. 

ARPER  *  Subindicator  (argument  of  perigee)  angular 

distance  of  perigee  from  node  at  reference  time. 
For  modified  orbital  elements,  this  is  also 
the  position  of  the  satellite  in  the  ellipse 
at  reference  time  (mean  anomaly  at  epoch  is 
always  equal  to  zero  in  this  system) 

11111  ■  angular  distance  of  perigee  from  northbound 
node,  measured  in  the  direction  of  satellite 
travel  in  degrees  and  hundredths  of  degrees 
m  ■  1  for  plus,  if  perigee  moves  in  the  same 

direction  as  satellite  travel.  2  for  minus,  if 
perigee  moves  in  the  direction  opposite  to 
satellite  travel 

nnn  »  P  dw/dt,  average  decimal  fraction  of  a  degree 
wfiich  perigee  moves  per  period,  measured  in 
thousandths  of  a  degree 
X  *  Always  an  X 

PERI0D  »  Subindicator  for  perigee- to-perigee  period 
(anomalistic  period) 

ooooo  ■  Anomalistic  period,  Pfl,  in  minutes  and  thou¬ 
sandths  of  a  minute,  if  first  two  digits  are 
less  than  85,  it  should  be  understood  that  100 
should  be  added  in  order  to  arrive  at  the 
correct  period.  Should  the  period  be  greater 
than  185  minutes,  a  special  notation  will  be 
made  in  the  message  . 

PPPPP  *  P  dP  /dt=-2c"P2  change  in  anomalistic  period 
pe  revolution,  measured  as  a  decimal  fraction 
in  one  hundred  thousandths  of  a  minute 
ECCEN  ■  Subindicator  for  eccentricity 
88888  ■  Eccentricity,  measured  as  a  decimal  fraction 
in  one  hundred  thousandths 

PERRA  ■  Subindicator  for  radial  distance  of  satellite 
from  center  of  earth  to  perigee 
rrrrr  ■  q,  radial  distance  of  satellite  from  center  of 
earth  at  perigee,  measured  in  miles  and  tenths 
of  miles 

RAFRE  ■  Subindicator  for  radio  frequencies  currently 
being  transmitted  from  satellite 
sssss  ■  Radio  frequency  in  megacycles  and  hundredths 
of  megacycles 

RTASC  *  Subindicator  for  right  ascension  of  the  node 
in  order  that  this  message  may  also  serve  the 
needs  of  those  who  prefer  traditional  orbital 
elements 
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ttttt  -  fl  ,  Degrees  and  hundredths  of  degrees 

Since  all  of  the  above  quantities  have  subindicators ,  a  message 
need  only  include  those  quantities  which  have  changed  since  the  last 
reference  time. 


SAT0R  code  output  is  computed  from  the  internal  elements  for  the 
reference  time,  the  time  of  perigee  passage,  T,  closest  to  epoch; 

M  ■  L  -0)  -  Cl  if  cos  i>0 

o  o  o  o 

M  *»  L  -03  +n  if  cos  i<0 

o  o  o  o 

T  -  t  -  -  M  /n 
o  o  o 

Compute  west  longitude  of  the  ascending  node  (T) ,  and 
argument  of  perigee ,  03  (T) ,  at  T:  fi 

9  (T)  (deg)  *  9  (at  start  of  epoch  year)  +  9  -  360°) 

gr  >  .  8^ 

D  +  9  f  where  9  is  the  rotation  rate  of  the  earth  in  deg/ 

solar  day,  D  and  f  are  T  in  days  and  fractions  of  a  day 

respectively  from  start  of  epoch  year. 

n<T>  -  -  'o' 

\  <T)  "egr(T)  "n(T) 

®  <T>  “  %  +  df  (T  ‘ 

The  "prime  sweep  interval"  is  the  mean  solar  time  interval 
between  the  passage  of  the  orbit  plane  over  a  geographic  location. 

Specif icially,  it  is  desirable  to  consider  the  recurrences  of  the  passage 
of  the  ascending  node  over  the  Greenwich  meridian.  Then  it  is  seen  that 
the  difference  of  the  mean  solar  time  of  this  event  from  one  day  to  the 
next, At,  is  due  to  two  causes:  the  revolution  of  the  Earth  about  the 
Sun  and  the  regression  of  the  nodes.  In  minutes  this  difference  is 


|(tji 


-  1 >  1440 


where  Sand  are  in  degrees  per  day. 
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APPENDIX  A 


SPADATS  ORBIT  REPRESENTATION 


The  following  formulation  is  shortly  to  be  used  by  SPADATS 
to  compute  ephemerides .  It  has  undergone  minor  changes  since  the 
latest  documentation  (SPADATS  Programming  Document,  Aeronutronic  Publi¬ 
cation  U-1691,  Revised  1  October  1962),  and  will  probably  change  further 
in  the  future. 

The  ephemeris  calculation  in  SPADATS  computes  a  mean  geocentric 
position  r  and  velocity  r.  at  some  time  t  from  the  epoch  orbit  elements 
a  N  ,  a  ,  t^,  and  Lq,  and  the  drag  parameter  co>  at  some  epoch  time  tQ. 
o  '  o  ° 

The  perturbations  caused  by  the  Earth's  atmosphere  and  by  the  asphericity 
of  the  Earth  (terms  due  to  the  second  and  third  zonal  harmonics)  are  in¬ 
cluded  in  this  calculation. 

The  constants  used  in  this  model  include  mean  equatorial  radius 
of  the  Earth: 


a  ■  6378165  meters 
e 

Second  Zonal  Harmonic  Coefficient 

J2  -  1.08245  x  10‘3 
Third  Zonal  Harmonic  Coefficient 

J3  »  -  2,50  x  10"6 

A.l  PRELIMINARY  COMPUTATIONS 

The  following  quantities,  required  in  later  calculations,  are 
first  determined. 
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(a)  Compute  the  semi-latus  rectum,  po,  the  eccentricity,  e  ,  the  mean 

semi-major  axis,  aQ,  the  perigee  distance  the  mean  value  of  the 

mean  motion  n  : 

o 


ao  (1  ‘  V 


(A-l) 


•fjr 


3/2 


eo2  (1  -  f  sin2  iQ) 


(b)  Compute  the  orientation  angles  iQ  and  : 


h 

-i  2 

i  *  cos”  ~~  ;  0<  i 
°  ° 


(used  for  output  only) 


n 

o 


the  quandrant  is  determined  from  the  signs 
of  the  numerator  (sine)  and  denominator 
(cosine) . 
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r 


a 


xN 


*yN 


.  CO.  <0  -  .  „  .in 

o  o 

*xN  Bin  “e  +  *yS  CO* 
o  o 


a 

—  ain  i 
P« 


o 


Tha  ■j*  term  in  a^  represents  the  leading  long-period  per¬ 
turbation  caused  by  the  Earth's  third  zonal  harmonic. 


(a)  Compute  the  mean  longitude  at  time  t,  as  perturbed  by  the  Earth's 
second  zonal  harmonic  (  Ait  term),  the  Earth's  third  zonal  harmonic 
(L^  term),  and  atmospheric  drag  (c  and  d  terms): 

(A-8)  L  •  L  +  n  (1+ A7T)  (t-t  )  +  L,  +  n  c  (t-t  )2  +  n  d  (t-t  )3 
OO  O  J  o  o  o  o 

where 

3  /  a» 

(A-9)  Att  -  |  J2  (  ^ 

(A-10)  L}..i  i  (j)  .1.  i, 

A. 3  CALCULATION  OF  GEOCENTRIC  POSITION  AND  VELOCITY  VECTORS 
(a)  Calculate  the  mean  argument  of  latitude  U  at  time  t: 

U  *  L  -  ft  if  i  <  90°  (direct  motion) 

U*L+  fti£i>  90°  (retrograde  motion) 


(2  -  |  sin2 


i  - 
o 


I  COS  i  I  ) 
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(b)  Solve  the  following  modified  form  of  Kepler'*  equation  for  the 
quantity  (E  +u> )  by  iteration,  uaing  a  first  guess  of  0  (mod  2^) 

(E  +co)  ■  U  +  a^  ain  (E  +ui )  -  a^  co*  (E  +to) 

(c)  Compute  the  mean  geocentric  position  r  and  velocity  r  at  time  t 
by  means  of  the  following  sequence  of  equations: 


e  cos  E  ■  a^  cos  (E  +u)  )  +  a^  sin  (E  +u> ) 


e  sin  E  ■  sin  (E  +a>)  -  ayN  cos  (E  +co) 


r  ■  a  (1  -  e  cos  E) 


i ,  iuJL 


e  sin  E 


•/Ma  J.  7 

rv  -  v  1  -  e 


C08  U  ■  — 


cos  (E  +C0  )  -  a^  +  ayN 


e  sin  E 


1  +  >  1 


sin  u  -  - 


sin  (E  +co  )  -  a  M  -  a, 


e  sin  E 


yN  “xN  j  +  ry 7 

V 


U  ■  N  cos  u  +  M  sin  u 


V  ■  •  N  sin  u  +  M  cos  u 


r  ■  rU 

r  -  *U  +  rv  V 
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APPENDIX  B 


DEFINITION  OF  OUTPUT  ELEMENTS 


This  appendix  defines  the  relationships  between  the  internal 
elements  used  by  the  SPADATS  Model  and  the  various  output  formats.  It 
also  presents  some  of  the  formulation  used  to  update  elements  after  a 
differential  correction. 


UPDATE  OF  ELEMENTS 


The  quantities  corrected  are  the  mean  argument  of  latitude,  U  , 
the  drag  coefficient,  c",  the  mean  value  of  the  mean  motion  at  the  epocfi, 
n  ,  and  the  parameters,  a^,  a  N,  Hand  1,  defined  in  the  text.  From 
the  new,  corrected  values  of  these  quantities  are  derived  the  period  decay 


rate. 


(B-l) 


c  ■ 

o 


-7? 

360 


(days/rev  ) 


and  the  mean  semi-n.ajor  axis  at  the  epoch, 

2/3 

V  ' 


(B-2) 


1  - 


J,a2 
2  e 


*  i  3,2 
(1  -  j  sin 


V 


J 


1  2/3 


1-e 


These  formulas  are  exact  inverses  of  equations  (A-l)  and  (A-2). 
Those  elements  experiencing  secular  variations  must  be  updated  to  the  new 
epoch,  t  ,  which  is  taken  at  the  crossing  of  ascending  node  nearest  the 
latest  observation.  As  in  Appendix  A,  the  rates  include 


(A-3) 


d  -  A(c")2 


An_  -  3n 
D  o 

3n_  -  3n 
D  o 


where  ■  0.072223521.  This  formula  is  analogous  to  one  found  by  Findley. 
The  constant  A*8  was  determined  by  experimentation*. 


/.  »  dn  j 

'A  5'  dt  “  ° 


3  Ve 

- r—  n_  cos  i 

4  O 


(A-6) 


do)  3  J2ae  r  .  ,  .  2  . 

dt  4  2  [  A  -  5  sin  i 

P« 


(A-9)  A7T  * 


.  2 
3  J2ae 
2  2 
P„ 


(2  -  ~  sin2  i  -  |  cos  i  |  ) 


The  most  complicated  update  formula  is  that  for  L 


■  K  +  n  (1  +  A7T)  (t  -  t  )  +  n  c"  (t  -  t  )2  +  n  d(t  -t  )3 
P  o  o  P  o'  o  vp  o  o  p  o' 


It  can  be  seen  that  the  periodic  term,  L_,  which  appeared  in  Eq.  A-8  is 
absent  because  thir  is  a  transformation  from  mean  elements  to  mean  elements 
at  a  new  epoch. 

The  argument  of  perigee  and  right  ascension  of  the  node  must 

be  updated 

40  ■  tan  * 
o 


*  For  details  see  Aeronutronic  Monthly  Status  Report  to  the  496L  SPO 
MPR-62-181,  October  1962,  p.  9. 


CD  ■  CO  +  (t  -  t) 

p  o  dt  p  o 


n  -  n  +  (t  - 1  > 

p  o  dt  p  o 


The  drag  parameters  cause  a  change  in  the  semi-major  axis 


(B-3) 


a  ■  a 

P  o 


During  the  calculation 
becomes 


+  2c"  (t  -  t  )  +  3d 
p  o 

of  the  new  elements, 


the  value  of  the  drag  coefficient 


(B-4)  cj  -  c"  +  3d  (tp  “  to) . 


FOUR-LINE  ORBITAL  ELEMENT  FORMAT  FOR  NORAD  SPADATS  USE 

The  orbital  element  standard  format  for  NORAD  SPADATS  use  has 
been  established  to  allow  Teletype  transmission  and  automatic  tape-to-card 
or  card-to-tape  conversion.  Line  0  is  used  only  for  descriptive  material 
of  primary  Interest  in  connection  with  the  issuance  of  bulletins.  Lines  1 
through  3  contain  all  the  numerical  values  of  the  orbital  elements  used 
in  computations  and  some  additional  values  of  interest  (the  anomalistic 
period,  the  perig< j  height,  and  the  rate  of  change  of  the  period)  which 
would  not  be  used  in  calculations.  Finally  a  check  sum  is  placed  at  the 
end  of  each  of  lines  1  through  3  to  allow  detection  of  transmission  errors. 

A  sample  message  as  it  would  appear  in  print  is  as  follows: 

0  005  312  58  BETA  002  US  11  62  04  27  314  04  29  98 

1  005  312  37781.00000000  01.35993597  -.168030-06  -.110000-10  120 

2  005  312  034.2490  275.4970  297.2870  .1862006  0133.99  000681  167 

3  005  312  156.2020  004.4050  -03.0140  -.612-04  -.140-3  016243  089 


LINE  0  SATELLITE  AND  BULLETIN  DESCRIPTION 


Col.  1 
Col.  2 
Cols.  3-6 


Col.  7 
Col.  8 


Cols.  9-11 
Col.  12 
Cols, 13-28 


Col.  29 
Cols. 30-33 


Line  Number  (0) 
Space 

Satellite  Number 
Col.  3 


Col.  4-6 
Space 

Source  of  Elements 


Element  Set  Number 
Space 

Astronomical  Name 

Cols .  13-14  Lest 

Col.  13  Space 

Cols.  16-24  Greek 

Col .  23  Space 

Cols.  26-28  Piece 

Space 

Country  of  Origin  and  Type 

Cols .  30-31  Code 

US 

SR 

FR 

UK 

CA 

JA 

Others 

Col.  32  Code 


Space  If  SPADATS  Number 

(Digit  4  through  9  If  SPASUR 
Number  used  prior  to  Issue  of 
SPADATS  Number) 

Satellite  Number 


Space  If  SPADATS 
9  If  SPASUR 


Lest  two  digits  of  year 
Space 

Greek  letter  (before  1963)  or  Numeral 
Space 

Piece  number  (before  1963)  or  Letters 


Country 
United  States 
Russia 
France 
England 
Canada 
Japan 

As  assigned 

Type 
Unknown 
Scientific, 
Weather 
Navigation 
Communications 
Manned 
As  assigned 
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LINE  1  SATELLITE  ORBITAL  INFORMATION 


Col.  1 
Cols.  2-12 
Cols.  13-26 


Col.  27 
Cols.  28-38 


Col.  39 
Cols.  40-50 


(B-5) 


Line  Number  (1) 

Same  as  L'.ne  0 

Epoch,  T,  Modified  Julian  Days 
(Julian  Day  minus  2,400,000.5) 

Cols.  13-17  Integral  part 

Col.  18  Decimal  point 

Cols.  19-26  Decimal  fraction 

Space 

Semi-major  Axis,  a,  Earth  equatorial  radii 
Cols.  28-29  Integral  part 

Col.  30  Decimal  point 

Cols.  31-38  Decimal  fraction 

Space 

First  Time  Derivative  of  Semi-major  axis, 
da/dt,  Earth  radii/day 

Sign 

Decimal  point 
Decimal  fraction 
Sign  of  exponent 
Exponent  of  10 


Col .  40 

Col .  41 

Cola.  42-47 
Col .  48 

Cols.  49-50 

^  t  a  c” 

dt  3  o 


Col.  51 
Cols.  52-62 


(B-6) 


pace 

acond_Time  Derivative  of  Semi-major  axis 
,-,J-  ,  Earth  radii/day/day 

52  Sign 

53  Decimal  point 

54.59  Decimal  fraction 

60  Sign  of  exponent 

61-62  Exponent  of  10 


*a/dt4 
Col. 
Col, 
Cols. 
Col. 
Cols . 


H  -  -  *  f 

dt* 


Col .  63 

Cols.  64-66 


Check  Sum  (Arithmetic  sum  of  digits  in  line, 
plus  one  for  each  minus  sign) 
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Col .  33 


Code 

0 

1 

2 

3 

4 

5-8 

9 


Type 

Silent 

Transmitting 
Rocket  body-silent 
Rocket  body-transmitting 

Metal  or  other  fragment 
As  assigned 

Unspecified 


Col.  34 
Cols .35-42 


Col.  43 
Cols. 44-46 
Col.  47 
Cols. 48-52 


Col •  53 
Cols. 54-56 


Space 

EpocU 

Cols.  35-36 
Col.  37 
Cols.  38-39 
Col .  40 

Cols.  41-42 
Space 

Bulletin  Number 


Last  tvo  digits  of  year 

Space 

Month 

Space 

Day 


"ace 

late  of  Bulletin  Issue 


Cols.  48-49 
Col .  50 

Cols .  51-52 
pace 

.ength  of  Bulletin 


Month 

Space 

Day 


Number  of  revolutions 
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LINE  2  SATELLITE  ORBITAL  INFORMATION 


Col.  1 
Cols.  2-12 
Cols. 13-20 


Col .  21 

Cols.  22-29 


Col.  30 
Cols.  31-38 


Col .  39 

Cols.  40-47 


Col.  48 
Cols.  49-55 


(B-7) 


Col.  56 
Cols.  57-62 


Col .  63 

Cols.  64-66 


Line  Number  (2) 

Seme  as  Card  0 
Inclination,  i.  Degrees 

Cols.  13-15  Integral  part 

Decimal  point 

Cols.  17-20  Decimal  fraction 

Space 

Argument  of  Perigee.oJ ,  Degrees 

Sis.  22-24  Integral  part 

_  ,  25  Decimal  point 

Coli.  26-29  °«l“l  t"ctl™ 

Right  Ascension  of  the  Ascending  Rode,  (1 

..  „  Integral  pert 

Decimal  point 
Decimal  fraction 


Cols.  31-33 
Col .  34 

Cols.  35-38 
Space 

Eccentricity, 

Col .  40 

Cols .  41-47 
Space 

Anomalistic  Period, 
Cols.  49-52 
Col.  53 
Cols.  54-55 

27T 

p  m  — — 
a  % 


(dimensionless) 

Decimal  point 
Decimal  fraction 


P  ,  Minutes 
“integral  part 
Decimal  point 
Decimal  fraction 


Height  of  Perigee  Above  Earth's  Equatorial 
Radius,  Hp,  Kilometers 

H  *  (a-ae-1)  a 
P  e 

Space  .. 

Check  Sum  (Described  in  Line  1) 


LINE  3  SATELLITE  ORBITAL  INFORMATION 


Col. 

Col8. 

Cols. 


Cols. 

Cols. 


Col. 

Cols. 


Col. 

Cols. 


(B-8) 


Col. 
Cols . 


1 

2-12 

13-20 


21 

22-29 


30 

31-38 


39 

40-47 


48 

49-55 


Degrees 

Integral  part 
Decimal  point 
Decimal  fraction 


Line  Number  (3) 

Same  as  Line  0 
Mean  Longitude, 

Cols.  13-15 
Col.  16 
Cols.  17-20 
Space 

First  Time  Derivative  of  Argument  of  Perigee, 
dco/dt,  Degrees  per  day 
Col.  22  Sign 

Cols.  23-24  Integral  part 

Col.  25  Decimal  point 

Cols.  26-29  Decimal  fraction 

Space 

First  Time  Derivative  of  Right  Ascension  of 
Ascending  Node,  dft/dt.  Degrees  per  day 
Col.  31  Sign 

Cols.  32-33  Integral  part 

Col.  34  Decimal  point 

Cols.  35-38  Decimal  fraction 

Space 

First  Time  Derivative  of  Eccentricity, 
de/dt,  per  day. 

Col.  40  Sign 

41  Decimal  point 

Decimal  fraction 
Sign  of  exponent 

Cols.  46-47  Exponent  of  10 


Col. 

Cols.  42-44 
Col.  45 


de 

dt 


5  (l-e) 


Space 

First  Time  Derivative  of  Anomalistic  Period, 
dP  /dt,  Minutes  per  day 
aCol.  49  Sign 

Col .  50  Decimal  point 

Cols. 51-55  Decimal  fraction  in  floating 

point  notation 
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dP 
_ a 


(B-9) 

—77  -  -  2c”  P 
dt  a 

Col. 

56 

Space 

Cols. 

57-62 

N  ,  Epoch  Revolution* 

Col. 

63 

Space 

Cols. 

64-66 

Check  Sura  (Described  in  Line  1) 

SEVEN' 

-CARD  ELEMENT 

FORMAT 

The  elements  used  In  the  formulation  described  in  Appendix  A 
are  introduced  into  the  computer  (and  produced  by  the  computer)  in  the 
form  of  seven  cards  shown  in  Figures  3-9.  Card  5  is  virtually  not  used 
by  the  system.  Card  7  contains  bookkeeping  information  concerning  the 
beginning  and  expiration  dates  of  bulletins  and  look  angles.  The  nodal 
period,  P  ,  on  card  6  is  calculated  by 


(B-10) 

P 

N  n 

1  -  !  J2  (; 

*eV  2  1 

(4-5  sin  i) 

0 

L  v 

«/ 

*  Launch  is  assumed  to  occur  in  Revolution  0.  Revolution  1  starts  at  the 
first  crossing  of  the  ascending  node.  The  revolution  count  is  not  always 
accurate  because  in  some  cases  revolutions  have  been  lost  due  to  lack  of 
sufficient  observations.  If  the  exact  time  of  launch  is  unknown,  an 
arbitrary  Revolution  1  may  be  defined. 
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Not  used 

A 
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- 
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E  ■  Nodal  Elements 


FIGURE  3  ELEMENT  CARD  1 
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FIGURE  4  ELEMENT  CARD  2 
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FIGURE  7  ELEMENT  CARD  * 
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FIGURE  8  ELEMENT  CARD  6 
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FIGURE  9  ELEMENT  CARD  7 
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RMS,  in  format  XXXXX.XX;  decimal 
punched  in  column  56 

Number  of  observations  used  in  ob¬ 
taining  RMS  IST0P 

Blank  or  0  -  correct  the  inclination 
element 

1  ■  do  not  correct  the  inclination 

2  ■  do  not  correct  the  drag  parameter 
4  ■  correct  time  equation  only 

Not  used 
Card  type 

E  -  Nodal  Elements 


APPENDIX  C 


DERIVATION  OF  PERTURBATIVE  TERMS 


For  the  actual  derivation  of  the  terns  in  Appendix  B  accounting 
for  the  asphericity  of  the  Earth's  potential,  the  reader  is  referred  to  the 
work  of  Kozai  (1959)*  and  of  Brouwer  (1959) .  The  intention  here  is  to  show 
that  the  SPADAT  System  formulas  agree  with  those  papers  except  in  notation. 

The  cited  papers  supply  some  of  the  forms  of  the  mathematical 
model,  the  constants  of  which  are  determined  by  fitting  the  observations. 
The  remaining  forms  are  designed  to  describe  the  effect  of  atmospheric  drag 
on  the  satellite.  The  mean  motion  is  described  as  a  quadratic,  and  the 
mean  longitude  is  represented  by  a  cubic.  The  mean  semi-major  axis  is 
related  to  the  mean  value  of  the  mean  motion  by  a  modified  form  of  Kepler's 
Harmonic  Law.  The  perigee  distance  is  assumed  unaffected  by  drag. 

The  formulas  of  the  Keplerian  ellipse  will  not  be  derived  here, 
even  though  their  formulation  in  terms  of  the  low-eccentricity  parameters 
may  be  unfamiliar  to  the  reader.  Thus,  this  appendix  will  be  restricted 
to  those  formulas  containing  the  perturbative  coupling  factors:  J_,  J. ,  c" 
and  d.  1  J 


KEPLER'S  HARMONIC  LAW 


Kepler's  Harmonic  Law  is  adopted  in  the  form  of  Kozai's  (1959) 
Equation  14,  which  in  the  notation  of  Appendix  A  is 


2  _ ' 

(1-  |  sin2  i  )  J  U? 
2  0*1  o 


★References  are  given  at  the  end  of  this  appendix. 
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It  must  be  remembered  that  the  subscripts  refer  to  (mean)  values  at  tQ. 


Kozai’s  a-  (unperturbed  semi-major  axis)  is  not  equivalent  to  aQ  (mean 
at  epoch),  above.  Equation  (A-l)  of  Appendix  A  follows,  if  terms  of 
order  J?  i 


are  omitted;  Equation  (B-2)  is  the  exact  inverse. 


The  translation  of  notation  will  be  facilitated  by  reference 
to  Brouwer's  (1959)  table  on  p.  396  and  by  comparing  the  potential  given 
on  that  page,  which,  in  SPADATS  notation,  is 


"•  V  [‘*£2  v£>J  v*i"s>] 


with  Kozai's  potential  on  p.  367: 


U- ~  ^1+“f  -  «in26>+  “f(f  *in2  6  -|)sin 


6  +. 


where  sin  ^ 


Thus, 


V  'J3  ae 


OTHER  SECOND- HARMONIC  TERMS 


The  nodal  regression  rate  and  movement  of  Perigee,  Equations 
(A-5  and  A-6)  are  derived  by  many  authors.  Kozai  shows  them  on  p.  372 
as  multipliers  of  elapsed  time.  Equation  (A-9)  is  the  sum  of  these 
rates,  with  the  unaer  standing  that  0  is  subtracted  in  retrograde  orbits 
and  that  nQ  is  separated  out  of  Att  (See  Equation  (A-9)  above). 


THIRD- HARMONIC  TERMS 


The  equivalent 
Here,  Kozai's  subscripted 
differences  between  e  and 
harmonics,  not  drag  as  in 


of  Equation  (A-7)  is  given  by  Kozai  on  p.  377. 
to  has  the  same  meaning  as  herein,  but  the 
e  are  periodic  terms  due  to  the  geopotential 
Equation  (A-7)  above.  Kozai  states 
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where 


"e  coscd  -  e  cos  to 
o 


—  3  A3 

e  sin  CD  *  e  sin  to+  7  ~ r  sin  i 
o  4  aA^ 


0)-  t  +  CD 
O 


The  perturbative  term  Is 


3  A3  .  .  1  J3  .  .  ®e 

7  —7  sin  i  ■  *  -  ~  sin  i  — 

4  aA2  2  J2  a 


2  This  is  equivalent  to  the  term  in  Equation  (A-7)  only  to 

order  e  .  Combination  of  Brouwer's  long-period  terms  in  e  and  g  -CD 

on  p.  394  will  also  give  Equation  (A-7)  and  show  that  the  effect  on 

a  „  is  of  order  e  . 
xN 


The  third-harmonic  term  in  the  mean  longitude,  L-,  Equation 
(A-10)  can  be  obtained  by  adding  the  three  formulas  of  Brouwer,  p.  391. 


L3  -  A3  i'  +  A3g'  +  A3h' 


where  the  lower  sign  corresponds  to  retrograde  orbits.  Jhus  in  Brouwer's 
notation 


r  1  .3  „  „  /  L'\  2  J  sin  I"  .  '  G"\ 
L3  “  4  7 2  c  8  \G"i  |  e"  1  "  [  L'J 

+  ^[hd_(H\2lp 

sin  I"  G"  \  G  /  f 


Translating  this,  but  not  combining,  gives 

L-  -  -  r  7“  ~  cos  (1-e  j 
3  2  J-  o  o 

Z  o 


1 

’o 


|cos  i  I  -  cos^i 


{sin  i 


1  -  (l"e0  ) 


By  expanding  and  omitting  the  eQ  terms,  one  obtains 


The  form  of  the  drag  perturbations  is  mainly  empirical.  The 
determination  of  the  observed  rate  of  change  of  the  orbit-period  from  the 
observational  data  appears  at  present  to  be  a  most  efficient  method  of 
predicting  future  behavior  of  the  satellite.  A  polynomial  expansion  of 
the  time  equation  such  as  used  for  many  years  by  SPACE  TRACK  (Findley  1962) 

VS)  WH  (»-»0)^0(N-»0)2«0(!)-H0)3 

o  o 


relates  the  nodal  crossing  time  tN  to  the  number  of  revolutions  (N-N  ) 
with  the  coefficients:  nodal  period,  rate  of  change  of  peiod  and  accel¬ 
erative  change  in  the  period,  P  ,  c,d,  respectively.  The  correction  of 
the  timing  through  the  expression  of  the  mean  longitude  has  been  more 
amenable  to  the  SPADATS  model.  It  is  necessary,  however,  to  establish 
the  relationship  between  the  foregoing  equation  and  Equation  (A-8) . 

In  all  that  follows,  the  J2  and  terms  will  be  omitted.  The 
coupling  of  these  terms  with  drag  is  negligible  for  most  satellites; 
where  it  is  not  negligible,  the  present  drag  theory  is  inadequate,  anyhow. 
This  inadequacy  is  confined  to  the  last  few  days  of  a  satellite's  life. 
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Derivatives  may  be  taken  of  the  old  SPACE  T*ACX  time  equation: 

i-  $  r  p» 

o  » 

Since 

.P  -1 
dt  N 

VPN  +2co<N'No)+3do(N-No)2 

o 

Similarly, 


ir  fvMo<"-"cp 


or,  approximately, 

2co"pn  3T 

From  the  mean  longitude  equation,  (A-8) , 


dL 


(c-D 


dfc  n-nQ  i  l+2c'*(t-to)+3d(t-t0) 


(C-2)  [lc"<«(t-t 0>  ] 

dt 


(C-3)  S-S 


f?  '6d  "■> 


The  anomalistic  period 


(B-7)  Pfi-  f" 

o 
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Is  equal  to  the  nodal  period  when  J 
Thus, 


2c 


p  fs 
N  dt 


dP 

7T 


2  ia  ignored. 

(2  7T  )2  dn  , 

2  dt 
n 


or 


(See  Equation  B-10.) 

-  SW?  2c"  „ 

3  o 

n 


The  factor  of  1440  between  the  last  equation  and  Equation  (B-l)  is  due 
to  the  fact  that  the  tine  is  measured  in  days  in  the  old  SPACE  TRACK 
equation  and  in  minutes  in  the  SPADATS  model. 

Equation  (A-3)  is  designed  to  permit  the  cubic  representation 
with  the  present  7-element  correction.  It  is  based  on  the  work  of 
Findley  and  must  be  accepted  as  purely  empirical.  The  relevant  exper¬ 
imentation  is  reported  in  Aeronutronic  Monthly  Status  Report  to  496L 
SPO,  MPR-62-181,  October,  1962,  p.  9.  Equation  (A-3)  may  be  replaced 
by  an  8-element  correction  soon.  The  formulas  of  Appendix  B  are 
consistent  with  either  source  of  d.  The  update  formula  (B-3)  is, 
again,  a  result  of  Kepler's  Harmonic  Law  in  the  form 

(C-4)  n-k  •JjT  a~3/2 
e 

and  Equation  (C-l) .  Formula  (B-4)  results  from  Equation  (C-2): 


(£)t  -no 

2c',+6d(tp-to> 

- 

(B-4)  cj;-c"+3d(tp-to) 

At  the  old  epoch  (or  the  new  epoch  after  the  update),  t*t 
(C-5)  ff  -  2c"  % 

From  E-iustion  (C-4), 


da 

2  a  dn 

2 

5  /  2  . 

a  dn 

dt  “ 

3  n  dt 

3 

k  dt 

e 
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The  first  form  gives 


(  da  t 

<»*»  \3t) 


-  I  ‘t  (2c"  "o>  -  -  3  *o  c" 


Differentiating  the  second  form  again,  results  In 


d2a 

dt2 


2  _a 

3 


d2n  5  a3^2  da  dn 

dt2  '  3  k  ^  dt  dt 


With 


the  help  of  Equations  (C-3),  <C-4)  and  <C-5>  ,  this  becomes,  at  tQ, 


(Si 


2  ®o 


f  -2  (6d  n  )  -  -r— 

3  n  o  3n 


-  -  A  a  d  +  a  (c*02 
o  7  o 


(' 1  v’) 


(2c”  n  ) 
o 


which  Is  Formula  (B-6) . 

From  the  assumption  that  the  perigee  height  does  not  change 


dt  d 


a  (1  -  e) 


-  0 


de  1-e  da 
dt  a  dt 


With  the  aid  of  Equation  (B-5) ,  this  becomes 


<»-»>  (fc)  -^7  (- 1  *. c") 

\  '  o  ' 


-  ?  (1  ‘  eo>  c" 
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